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SUMMARY 


The  problem  of  the  ‘build--up  of  lift  on  two—  and  three-dimensional 
wings  flying  at  high  speeds  is  discussed  as  a  boundary— value  problem  for 
the  classical  wave  equation.  Kirchhoff  *s  formula  is  applied  so  that  the 
analysis  is  reduced,  just  as  in  the  steady  state,  to  an  investigation  of 
sources’  and  doublets.  Some  simple  applications  of  this  method  are  con¬ 
sidered,  including  the  determination  of  the  starting  lift  of  a  three- 
dimensional  wing  and  the  potential  functions  for  some  types  of  -unsteady 
vortex  motion. 

INTEODUCTIOK 


The  usual  idealizations  Introduced  in  the  development  of  linearized 
aerodynamic  theory  describe  a  frictionless,  perfectly  elastic,  model 


fluid, 

fluid 

terns 

,  As  is  well  known,  the  effect  of  small  disturbances  in  such 
can  be  analyzed  by  means  of  the  familiar  wave  equation  which, 
of  the  pertiu"bation  velocity  potential  q),  can  be  written 

a 

in 

^xx  ^yy  ^zz 

Vff 

(1) 

where 

ao 

speed  of  pressure  propagation 

t* 

time 

spatial  coordinates 

3*  Equation  (l)  applies’  to  a  flow  field  which  is  stationary  at  large  dis¬ 

tances  from  the  disturbance  region;  furthemore,  the  axial  system  is 
stationary  relative  to  the  fluid  infinitely  distant.so  that,  if  a  moving 
%  wing  is  being  analyzed,  the  wing  moves  with  respect  to  the  x,y,z  axes. 


It  is  possible  to  use  a  transformation  which  will  bring  the  wing  to 
rest  with  respect  to  a  new  axial  system^  but  this  process  leads  generally 
to  a  different  equation  for  the  perturbation  potential.  The  case  of  a 
wing  moving  at  subsonic  speed  constitutes  an  exception^  for  here  the 
Lorentz  transformation^  familiar  in  modem  physics^  will  fix  the  origin 
of  the  axial  system  on  the  wing  while  leaving  the  wave  equation  Invariant, 
Such  is  not  the  case  when  the  wing  travels  at  supersonic  speeds,  so  that, 
if  it  is  desired  to  fix  the  axial  system  to  the  wing,  the  wave  equation 
must  be  abandoned.  This  latter  alternative  has  been  the  more  extensively 
studied  in  aerodynamics,  but,  because  of  the  large  body  of  knowledge  con¬ 
cerning  equation  (l)  which  is  available  in  mathematical  physics,  the 
present  report  is  based  upon  the  wave-equation  approach  to  problems  of 
unsteady  motion. 

The  solution  to  equation  (1)  known  as  Kirchhoff*s  formula  (see  ref¬ 
erence  1)  is  found  to  be  of  considerable  use  in  unsteady-motion  problems 
involving  thin  wings  with  supersonic  edges.  The  problem  is  reduced  to 
one  of  simming  elementary  solutions,  analogous  to  sources  and  doublets 
in  steady  flow,  over  a  region  determined  by  the  position  of  the  wing  as 
well  as  its  traversed  path.  The  theoretical  development  leads  naturally 
to  the  concepts  (defined  later)  of,  first,  inverse  sound  waves,  which 
have  a  counterpart  in  the  Mach  forecones  used  in  steady— state  wing  theory; 
second,  acoustic  plan  forma;  and,  third,  homogeneous  flow,  which  reduces 
in  part  to  the  familiar  conical  flow  as  the  wing  approaches  a  steady 
supersonic  velocity. 

In  connection  with  boundary— value  problems  involving  noninteracting 
surfaces,  a  theorem  will  be  given  which  makes  it  possible  to  build  up  a 
three-dimensional  supersonic— edged  wing  by  superposition  of  elements  of 
two-dimensional  wings.  This  result  is  an  extension  to  unsteady  flow  of 
a  theorem  given  previously  by  Lagerstrom  and  Van  Dyke  (reference  2).  By 
means  of  this  theorem,  problems  involving  super sonlc-edged  wings  per¬ 
forming  any  prescribed  maneuver  can  be  solved,  provided  only  that  the 
behavior  of  a  two-dimensional  wing  executing  the  same  maneuver  is  known. 
Thus,  the  relatively  large  amount  of  material  available  for  two-dimensional 
unsteady  motion  at  supersonic  speeds  can  be  carried  over  directly  to  the 
supersonic-edged  wing. 

In  the  last  part  of  the  report,  some  simple  applications  of  the  gen¬ 
eral  methods  are  considered,  including  the  determination  of  the  starting 
lift  of  a  three-dimensional  wing  and  the  potential  functions  for  some 
types  of  -unsteady  vortex  motion. 


UST  OF  IMPOETANT  SYMBOLS 


aQ  free-stream  speed  of  soimd 

Cq  vlng  root  chord,  or  maximum  chord  of  a  wing 
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c 

L 

Mo 

q 


ac. 


wing  chord 

section  lift-curve  slope 
lift  coefficient 

coefficient  of  rolling  moment  due  to  yawing 
lift 

free-stream  Mach  number 


loading  coefficient  (pressure  on  the  lower  surface  minus  pressure 
on  the  upper  surface  divided  by  dynamic  pressure) 


a(rt/2Voy 


q 

t» 

t 

w 


( I 


free— stream  dynamic  pressure 
time 
aot’ 

free-stream  velocity 

vertical  component  of  perturbation  velocity 
XyjyZ  Cartesian  coordinates 
a  angle  of  attack 

p  y  |i  ^Mo^i 

\  slope  of  stream  surface  in  free-stream  direction 

p  free-stream  density 


cp  perturbation  velocity  potential 

A  discontinuity  in  the  quantity  in  question  across  the  z  =  0  plane 

Subscript 


u  value  on  upper  surface  (z  =  0  plane)  of  a  wing 
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THEORETICAL  DEVELOPMENTS 


Basic  Equation  and  Solutions 


The  partial  differential  equation  which  foms  the  basis  of  the  pres¬ 
ent  study  has  been  written  as  equation  (l).  It  is  possible  to  put  this 
equation  in  a  more  convenient  form  by  Introducing  the  notation 


t  =  aot» 


so  that  the  dimension  of  t  is  length  Just  as  are  the  dimensions  of  the 
geometric  variables  x,  y,  and  z.  Equation  (2),  together  with  equation 
(1),  yields  the  canonical  form  of  the  wave  equation 


‘Pxx  +  'Pyy  +  'f’zz  =  'Ptt  ^3) 

and  it  is  this  form  which  will  be  considered. 

The  first  task  is  to  study  the  relation  between  the  motion  of  the 
wing  and  the  coordinate  system.  As  has  already  been  mentioned^  equa¬ 
tion  (3)  is  valid  for  a  flow  field  produced  by  a  wing  moving  relative  to 
a  fixed  coordinate  system.  It  is  pertinent  to  consider  the  possibility 
of  finding  a  transformation  which  will  (l)  fix  the  origin  of  the  axial 
system  on  the  wing  and,  at  the  same  time,  (2)  retain  the  wave  equation 
as  the  governing  equation  of  the  flow.  Certainly  the  first  of  these 
requirements  is  simple  to  fulfill  if  the  second  is  neglected.  However, 
because  of  the  simplicity  of  the  wave  equation  and,  what  is  more  impor¬ 
tant,  because  of  the  great  amount  of  developmental  study  that  has  been 
expended  on  it,  the  second  reqiilrement  is  not  without  Justification. 

The  following  transformation  (known  in  relativity  theory  as  a  Lorentz 
transformation,  or  as  a  hyperbolic  rotation) 

x-Mpt 

^  ■  /l-Mo^ 

n  =  y 

C  =  z 

t-MpX 


>  W 
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where  !!„  =  — 
a, 


<1,  will  satisfy  "both  conditions  (l)  and  (2),  For  example. 


suppose  that  the  wing  is  moving  along  the  x  axis  with  velocity  Yq. 
Application  of  eqiiation  (4)  to  equation  (3)  makes  the  origin  of.  the  new 
axial  system  also  travel  along  the  x  axis  with  a  velocity  Yq.  This  is 
seen  to  he  so  since  |  is  always  zero  when  x  =  Vot'  =  Mot.  Hence,  the 
I  axis  is  "fixed"  on  the  wing.  As  to  the  second  condition,  a  straight¬ 
forward  exercise  in  partial  differentiation  yields 


'Ptt  ^ 

so  that^  in  going  from  to  space^  the  wave  equation 

remains  invariant;  consequently^  toth  the  requirements  mentioned  have 
been  fulfilled. 


It  is  Instructive  to  consider  briefly  the  consequences  of  applying 
the  Lorentz  transfoimiation.  Although  the  wave  equation  remains  invariant, 
such  physical  quantities  as  length  and  pressure  do  not.  For  example,  a 
wing  with  a  chord  Cq  in  the  x,y,z  space  has,  according  to  equation  (4), 

a  chord  Cq/v^  1— in  the  space.  Furthermore,  the  loading 

coefficient  which,  on  the  basis  of  linearized  theory,  is  given  in  the 
x,y,z,t  space  by 


Ap  _  4  39 

q  YqMo  St 

becomes  for  the  |,  t},  t  space 

Ap _ h  /  1  ^(p  S<p\ 

^  ~  VqV^I-  Mo  2  St  J 


(6) 


(7) 


If  the  wing  motion  is  steady  and  there  are  no  transient  effects,  equa¬ 
tions  (5)  and  (7)  are  independent  of  time  and  together  with  the  resiilting 
length  transformations  become 

'Pll  +  9^^  =  0 

Ap 


> 


(8) 
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These  are  immediately  recognized  to  1©,  respectively,  Laplace’s  equation 
for  incompressihle  flow  and  the  familiar  Prandtl-Glauert  compressibility 
corrections.! 

The  preceding  discussion  has  an  important  qualification,  however,  in 
the  fact  that  the  velocity  of  the  moving  axial  system  cannot  exceed  the 
speed  of  sound,  A  glance  at  equation  (4)  serves  to  verify  this  statement 
since  that  equation  shows  Mq  must  be  less  than  1  in  order  that  i  and 
T  b.e  real  for  real  x  and  t.  In  fact,  it  has  been  shown  that  there  is 
no  transformation  which  will  fix  the  moving  axial  system  in  a  wing  travel¬ 
ing  at  a  uniform  supersonic  speed  away  from  the  original  fixed  axes  and 
still  keep  the  wave  equation  invariant.  Therefore,  for  analyzing  a  wing 
in  supersonic  flight,  it  is  necessary  to  abandon  one  of  the  two  proposed 
requirements:  either  the  axial  system  cannot  be  fixed  in  the  wing,  or 
the  field  equation  must  be  modified.  The  latter  of  these  two  alternatives 
has  been  studied  by  several  authors  (e.g.,  references  3,  h,  5,  6)  but  it 
is  the  fomer  which  will  be  considered  in  the  present  analysis.  Further, 
since  the  axes  cannot  be  made  to  travel  as  fast  as  the  wing,  they  will 
not  be  made  to  move  at  all  and  equation  (3)  will  be  adopted  throioghout  as 
,  the  basic  equation. 

Having  decided  upon  the  form  of  the  partial  differential  equation, 
the  boundary  conditions  must  next  be  established.  For  any  given  time 
these  conditions  are  similar  to  those  studied  in  steady— state  thin— 
airfoil  problems;  namely,  either  that  the  given  slope  of  the  wing  surface 
is  proportional  to  the  vertical  induced  velocity  cp^  over  the  region 
occupied  by  the  wing  in  the  z  =  0  plane, ^  or  that  the  prescribed  surface 
pressure  is  proportional  to  the  timewise  gradient  cp^  in  velocity  poten¬ 
tial  over  the  same  region.  The  addition  of  time  simply  means  that  this 
region  moves  about  in  the  z  =  0  plan©  in  conformity  with  the  known 
direction  and  velocity  of  the  wing. 

The  additional  condition  is  imposed  that  the  induced  velocities  fall 
to  zero  on  the  surface  which  is  formed  by  the  envelopes  of  the  spherical 
sound  waves  originating  from  the  surface  of  the  wing. 

The  solution  to  equation  (3),  subject  to  the  boundary  conditions 
Just  mentioned,  can  be  expressed  by  a  formula  which  may  be  regarded  as 

!The  equation  |=  x/^ T-Mq^  would  first  read  |  =  (x-MotQ )  / 

where  to  is  a  constant  representing  the  time  required  for  the  motion 
to  reach  its  steady  state.  However,  the  x  coordinate  can  always  be 
translated  to  any  fixed  position  without  affecting  any  of  the  equations 
for  potential,  loading,  etc.  Such  a  translation  is  assumed  to  have 
been  mad©  in  equation  (8) . 

^The  z  =  0  plane  is  assumed  to  be  the  "plane  of  the  wing";  that  is,  if 
the  angle  of  attack  were  zero  and  the  wing  had  no  thickness  it  would 
lie  entirely  in  the  z  =  0  plane. 
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requiring  either  the  evaluation  of  a  double  integral  or  the  solution  of 
a  double -integral  equation,  depending  upon  whether  a  boundary— Value 
problem  of  first  or  second  kind  is  considered.  This  solution  is  known 
as  Kirohhoff's  formula.  (See  reference  1.)  It  may  be  written  in  a  form 
convenient  for  aerodynamic  applications  as  follows: 


_  A  A  —  - j - J 


dS  (9) 


’a 


where  the  A  indicates  the  Jump  (value  on  the  upper  surface  minus  value 
on  the  lower  surface  when  applied  to  cp  or  of  the  function  in 

passing  through  the  2i=0  plane,  r  =  xhki ) ^+( y— yi ) ^+( z~zi ) ^  j 

^o  “  a/( ^““^1 )  y  71 ) ^+2 ^  ,  and  where  the  area  of  integration  Sa  will  he 
discussed  in  more  detail  later. 


The  terms  in  the  Integrand  of  equation  (9) 
introducing  the  following  notation: 


can  he  shortened  hj 


a 


<p(xi,yi,zi,t-^o)  =  ^  ^ 


(10a) 


Szi  r 


■  dr  d  9(xi,yi,0,t-r) 

=  A-1 

'i' 

_dza  dr  r  _ 

bz  br 

r 

(lOb) 


In  this  notation,  the  subscript  r  in  equation  (lOa)  means  that  r  is 
to  he  held  constant  in  the  differentiation,  and  the  prefix  A  obviates 
the  necessity  of  Indicating  that  the  functions  considered  are  to  he 
evaluated  for  Zi=0,  since  it  Indicates  that  the  difference  of  the  values 
of  the  function  across  the  Zi=0  plane  is  to  he  taken.  The  right-hand 
side  of  equation  (iOa)  can  he  recognized  as  a  term  representing  a  source 
located  in  the  Zi=0  plane,  and  the  right-hand  side  of  equation  (lOb)  Is 
seen  to  represent  a  doublet  located  in  and  with  axis  normal  to  the  zi=0 
plane.  The  brackets  [  ]  about  the  functions  in  equations  (lOa)  and 
(lOb)  have  a  special  meaning  which  is  defined  in  the  following  way:  if 
f  is  a  function  whose  value  at  a  fixed  point  P  depends  upon  the  coor¬ 
dinates  xi,yi,zi,t  of  a  moving  point  Q,  so  that 

f  =  f(xi,yi,zi,t) 


[f]  =  f(ii,yi,zi,t-r) 


(11) 


then 
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where  r  is  the  distance  from  P  to  Q,  As  an  example,  consider  the 
potential  cp  at  a  point  P  due  to  a  moving  source,  the  location  of  which 
at  any  time  is  Q.  Then  cp  satisfies  the  condition  just  mentioned  that 
it  depends  on  the  coordinates  Xi,yi,Zi,t  of  Q.  The  "brackets  [  ]  indi¬ 
cate  that  the  potential  [  cp  ]  depends  not  upon  the  source  strength  now  at 
"time*^  t,  hut  rather  upon  the  source  strength  that  existed  "time"  r 
ago.®  For  convenience,  [cp]  is  referred  to  as  the  retarded  value  of  cp. 

The  expression  for  a  doublet,  equation  (lOh),  is  usually  expanded 
as  follows; 


■9  ' 

^  A 

A  I”’ 

Bz  Sr 

. 

Sz  Si’oVi’c 

J  Sz  ro 

Sr  ^ 

Finally,  equation  (9)  becomes 


5t 


cp  = 


4it 


rn 


5Acp 
Szi  J 


fo  dz 


^A(p 

L 


} 


dS 


(12) 


(13) 


The  application  of  equation  (13)  awaits  only  a  discussion  of  the  area  Sa 
over  which  the  integration  is  to  he  made.  This  discussion  is  important 
enough,  however,  to  merit  consideration  in  some  detail  and  will  be  given 
in  the  following  section. 


The  Acoustic  Plan  Form 


Suppose  that  a  line  of  sources  is  placed  along  the  y  axis  and  that 
the  strength  of  these  sources  is  zero  for  t  <0.  At  t  =  0  they  are 
"turned  on"  and,  at  the  same  time,  start  moving  along  the  negative  Xi 
axis  with  the  velocity  Vq-  After  time  t'  has  passed,  the  source  line 
has  traveled  a  distance  l^t  as  shown  in  the  accompanying  sketch  (which 
is  drawn  for  the  case  Mq>1). 


3 

Quotes  are  used  around  the  word  time  since  the  dimension  of  t  is 
actually  length,  not  time.  It  is  convenient,  however,  to  refer  to  t 
as  "time,"  and,  since  the  actual  value  of  time  is  simply  t  multiplied 
by  the  constant  a^,  this  should  cause  no  confusion. 
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Suppose  next  that  there  are  tvo  sensing  elements,  or  detectors,  placed  at 
the  point  P(x,y)  located  somewhere  ahead  of  the  axis;  one  of  these 

detectors  is  responsive  to  light  and  the  other  to  sound,  Now,  the  light 
detector  at  any  given  time  will  show  the  sources  lying  in  a  straight  line 
just  as  they  would  appear  visually  at  each  particular  instant.  The  situa¬ 
tion  is  entirely  different,  however,  for  the  sound  detector.  First  it  is 
necessary  to  understand  the  nature  of  a  spherical  sound  wave.  Such  a 
wave  travels  outward  from  its  origin  at  a  velocity  ao>  so  that  in  the 
time  t*  it  has  traveled  a  distance  t.  Before  the  wave  reaches  a  point, 
the  point  is  completely  unaware^  of  its  existence  and,  further,  after  the 
wave  has  passed,  the  point  remains  subsequently  \maware  of  its  existence. 
Hence  the  only  points  disturbed  by  the  wave  are  those  momentarily  on  the 
spherical  surface  Itself ,  (in  this  connection  see  reference  1,  pp,  1  ff.) 


.4 

IMaware  is  used  in  the  sense  that  an  Instrument  will  record  no  change 
in  any  of  the  physical  properties  of  the  air  at  the  point  in  question. 


10 


MCA  TIT  2256 


The  sound  detector,  therefore,  can  only  "hear"  sources  which  are  so 
located  that  their  spherical  soimd  waves  are  Just,  at  the  given  instant 
of  time,  reaching  the  detector.  The  locus  of  all  the  points  which,  at  a 
time  t  ago,  emitted  sound  waves  that  are  just  now  reaching  the  point 
P(x,y)  is  itself  a  sphere  and  for  convenience  this  sphere  will  he 
referred  to  as  an  inverse  sound  wave.^  The  traces  of  these  inverse  sound 
waves  In  the  z  =  0  plane  are  drawn  in  the  sketch  as  concentric  circles 
about  the  point  P(x,y).  The  intersection  of  an  inverse  sound  wave  of 
radius  t— t  with  the  line  representing  the  position  of  the  sources  at  a 
time  t*— ago  gives  the  position  of  the  sources  which  are  Just  now 
signaling  their  presence  to  the  sound  detector  at  P,  For  example,  when 
the  source  line  started,  it  was  lying  along  the  yi  axis.  With  referenc 
to  the  present  time  t  this  was  aot’  removed.  Hence  the  intersection 
of  the  circle  about  P  of  radius  t  with  the  yi  axis  fixes  the  two 
points  A  and  A*,  the  spherical  sound  waves  of  which  are  now  reaching 
P,  A  continuation  of  this  process  yields,  for  the  locus  of  all  points 
from  which  waves  emanated  that  are  Just  now  touching  the  point  P,  the 
part  of  the  ellipse  shown  in  the  sketch. 

In  the  sense  that  the  light  detector,  because  of  the  very  large 
velocity  of  light,  is  "seeing"  a  straight  line  of  sources,  the  sound 
detector,  because  of  the  relatively  slow  velocity  of  sound,  is  "hearing" 
an  elliptic  line  of  sources.  Extending  this  concept  to  include  a  sheet 
of  sources  distributed  over  the  surface  of  the  wing,  one  can  refer  to 
the  outline  of  that  part  of  the  wing  which  generates  disturbances  which 
can  be  measured  by  the  light  detector  as  the  plan  form  (l.e.,  the  visual 
plan  form),  and  to  the  outline  of  that  part  of  the  wing  which  affects 
the  sound  detector  as  the  acoustic  plan  form.  In  a  mathematical  sense, 
the  acoustic  plan  form  is  the  area  S  over  which  the  integration  of 
equation  (I3)  is  to  be  made. 

The  equation  for  the  acoustic  plan  form  can  be  formulated  by  means 
of  the  two  equations 


+  (y-yi)^  +  =  (t-r)^ 

(14) 

f(yi,xi,T)  =  0 

(15) 

where 

x,y,z  coordinates  of  the  point  at  which  the  induced  effects  are  to  be 
measured 

^The  inverse  sound  wave  has  for  its  analogue  in  steady  lift ing-eurf ace 
theory  the  Mach  forecone.  In  that  theory  a  disturbance  outside  the 
Mach  forecone  cannot  affect  the  values  of  the  Induced  velocities  at 
the  point  where  they  are  being  measured.  Similarly,  in  the  present 
study,  a  source  located  outside  the  inverse  sound  wave  of  radius  t 
cannot  affect  the  values  of  any  measurement  made  at  the  point  P, 
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^ijyi  variable  points  of  the  sources 
t  "time"  now 

t— T  "time"  ago 


Equation  (l4)  Is  that  of  the  Inverse  sound  waves  and  equation  (15)  repre¬ 
sents  the  position  of  the  visual  plan  foim  at  a  time  t— t  past.  It  Is  nec¬ 
essary  to  Include  the  region  hehlnd  the  wing  covered  hy  the  vortex  wake 
as  part  of  the  visual  plan  form.  In  case  the  vortlclty  In  the  wake 
vanishes,  as  In  the  thickness  problem,  the  wake  may  still  he  considered 
as  part  of  the  visual  plan  form,  hut  the  strength  of  the  source-douhlet 
distribution  over  that  part  of  the  acoustic .plan  form  corresponding  to 
the  wake  will  vanish.  If  t^— z®<0,  the  acoustic  plan  form  does  not 
exist  since,  for  such  a  case,  no  source  on  the  Xi,yi  plane  has  had  time 
to  transmit  Its  effect  to  the  point  x,y,z.  On  the  other  hand.  If  the 
circle  (in  the  Xi,yi  plane)  given  hy  the  equation 

(x-x,)^  H-  (y-yi)^=  t^-^  = 


lies  entirely  within  the  area  occupied  hy  the  visual  plan  form  of  the  wing 
at  the  beginning  of  the  motion,  the  acoustic  plan  form  Is  just  this  circle 
Itself.  For  any  other  situation  the  acoustic  plan  form  Is  formed  In  part, 
or  In  whole,  hy  the  curve  found  hy  eliminating  t  from  equations  (li|-) 
and  (15)  and  In  part,  or  not  at  all,  hy  an  arc  of  the  circle  (x— Xi)^  + 
(y~yi)^  =  t®-2^.  One  of  the,  principal  advantages  of  the  acoustic— plan- 
form  concept  arises  In  the  study  of  problems  Involving  source  or  doublet 
distributions  having  constant  strength.  In  such  cases,  the  retarded 
values  of  the  potential  and  Its  gradient  appearing  In  equation  (13)  are 
constant  and  can  be  taken  outside  the  Integral  signs.  The  problem  la 
thereby  reduced  to  the  Integration  of  a  simple  geometric  variable  over 


A  few  examples  will  serve  to  fix  the  Idea  of  the  acoustic  plan  foim. 
Consider  first  a  two-dimensional,  unswept  wing  moving  at  a  constant  super¬ 
sonic  speed  In  the  negative  Xi  direction.  At  time  zero  the  leading  edge 
of  the  wing  was  along  the  yi  axis  and  now,  at  time  equal  to  t*,  the 
wing  has  moved  so  that  the  leading  edge  coincides  with  the  line  Xi  =  — 

Mot .  Choose  three  points  that  are  now  lying  on  the  wing .  Let  one  point 

have  its  x  coordinate  In  the  range  Cq— \/t  2-2;2  >  X>yt2-z2 

(where  Cq 


is  the  chord  of  the  wing),  the  second 

y t^— z^  ,  and  the  third  is  the  range 
these  points  by  Pi,P2,  and  Pg  (see 
acoustic  plan  forms  are,  respectively, 
part  ellipse,  and  a  complete  ellipse. 


in  the  interval 


y  t^-z 


>x> 


— yt^— z^  ^  X  ^  “MQ't*  Designating 
sketch).  It  can  be  shown  that  their 
a  complete  circle,  a  part  circle  and 
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The  points  P  are  at  the  centers  of  the  circles  and  at  focal  points  of 
the  ellipses.  Since ^  moreover,  the  circTilar  plan  form  about  Pi  receives 
’no  signals  from  sources  on  the  leading  or  trailing  edge,  conditions  at  Px 
are  consequently  completely  Independent  of  the  actual  (visual)  plan  form 
of  the  wing.  The  elliptical  plan  form  about  P3,  on  the  other  hand, 
depends  entirely  on  the  shape  of  the  leading  edge;  and  finally  the  mixed 
plan  form  about  P^  is  in  certain  regions  (the  circular  portion)  independ¬ 
ent  of  the  leading  edge,  and  in  other  regions  (the  elliptic  portion) 
entirely  dependent  upon  it.  Since  the  wing  is  traveling  at  supersonic 
speeds,  the  trailing  edge  and  vortex  wake  can  have  no  effect  on  the  meas¬ 
urements  taken  on  the  wing  and,  in  the  same  way,  a  point  ahead  of  the  wing 
leading  edge,  P^  in  the  sketch,  is  undisturbed. 

Next  consider  a  wing  moving  at  a  constant  subsonic  speed  in  the  nega¬ 
tive  Xi  direction.  As  before,  the  leading  edge  was  on  the  yi  axis  at 
t*  =  0  and  has  traveled  a  distance  Choose  now  three  points  Pi, 

Pg^  and  Pg  on  the  wing  and  unaffected  by  the  wing  tips.  The  acoustic 
plan  forms  for  these  points  are  combinations  of  circles  and  hyperbolas  as 
contrasted  with  the  circle— ellipse  combination  in  the  supersonic  case. 

Just  as  in  the  supersonic  case,  however,  there  is  a  certain  region 
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represented  ty  Pi  in  which  the  acoustic  plan  form  is  a  complete  circle 
and  is  independent  of  the  visual  shape  of  the  wing  (see  sketch). 


Xl 


Point  P2  is  surrounded  hy  a  plan  form  which  is  part  hyperbolic  and  part 
circular,  the  point  itself  being  the  center  of  the  circle  and  the  focus 
of  the  hyperbola.  Point  P3  is  a  limiting  value  of  P2;  it  lies  on  the 
leading  edge  of  the  wing  and  the  hyperbolic  sides  of  its  plan  form  have 
degenerated  into  straight  lines .  Finally,  P^  lies  ahead  of  the  wing; 
its  plan  form  is  still  a  combination  of  a  hyperbola  and  a  circle,  but  P4 
is  now  the  focal  point  lying  ahead  of  the  hyperbolic  branch  used. 

The  sketch  was  constructed  so  that  the  portion  of  the  visual  plan 
form  behind  the  trailing  edge  had  no  effect  on  the  potential  at  the  var¬ 
ious  points  Pi,  etc.  If  these  points  had  been  chosen  at  positions  where 
the  wake  could  signal  its  effect,  one  of  two  acoustic  configurations  would 
result.  First,  if  the  wing  is  symmetric  about  the  z  =  0  plane,  no  lift 
is  developed  and  the  vorticity  in  the  wake  is  zero  so  that  the  visual  plan 
form  need  not  include  the  wake,  but  effectively  ends  at  the  trailing  edge. 
In  this  case,  the  leading  edge  of  the  acoustic  plan  form  is  then  deter¬ 
mined  as  before,  while  its  modified  trailing  edge  may  be  made  up,  in  part. 
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of  circular  arcs  formed  by  the  primary  wave  and,  in  part,  by  an  arc  of  the 
hyperbola  formed  by  the  (acoustic)  intersection  of  the  straight  visual 
trailing  edge  with  the  primary  wave  (such  an  arc  being  identical  with  the 
leading  edge  of  the  acoustic  plan  form  but  displaced  backwards).  On  the 
other  hand,  if  the  wing  has  no  thickness  but  is  inclined  to  the  free 
stream,  it  develops  lift  and  the  vorticlty  in  the  wake  does  not  vanish; 
the  acoustic  plan  form  has  a  trailing  edge  made  up  entirely  of  an  arc  of 
the  primary  inverse  sound  wave.  The  space  between  this  arc  and  the  acous¬ 
tic  trace  of  the  visual  trailing  edge  is  covered  by  a  sheet  of  doublets, 
the  strength  of  which  is  determined  by  the  vortlcity  distribution  of  the 
vortex  wake. 

It  is  interesting  to  notice  the  conversion  of  terminology  which 
arises  in  the  analysis  of  unsteady  lift  problems.  In  the  study  of  steady- 
state  wings,  it  is  customary  (because  of  the  nature  of  the  governing 
partial  differential  equation)  to  speak  of  the  subsonic  problems  as  ellip¬ 
tic  and  the  supersonic  problems  as  hyperbolic.  Yet  the  acoustic  plan 
forms  Just  presented  involved  ellipses  for  the  supersonic  wing  and  hyper¬ 
bolas  for  the  subsonic  case. 

To  complete  the  remark,  it  can  be  observed  that  when  the  velocity  of 
the  wing  is  sonic  the  steady— state  partial  differential  equation  becomes 
parabolic  and,  in  this  case,  the  acoustic  plan  form  of  a  straight-edged 
wing  also  involves  parabolas.  When  the  leading  edge  is  linear  and  normal 
to  the  stream  direction,  the  eccentricity  of  the  conic  sections  bounding 
the  acoustic  plan  form  for  a  point  on  the  wing  is  equal  to  I/Mq.  From 
this  relation  it  is  apparent  that  for  Mq  less  than,  equal  to,  and  greater 
than  1  the  sections  are,  respectively,  hyperbolas,  parabolas,  ajid  ellipses. 
As  might  be  presumed,  the  value  of  the  eccentricity  satisfies  simple  sweep 
theory  so  that,  for  an  infinitely  long  straight  leading  edge,  the  eccentric¬ 
ity  of  the  acoustic  plan  form  is  I/Mq  cos  A  where  A  is  the  angle  of 
sweepback.  The  principal  axes  of  the  conic  sections  are  always  normal 
and  parallel  to  straight  leading  edges. 


Homogeneous  Boundary— Value  Problems 


Kirchhoff *s  solution  to  the  wave  equation  can  be  applied  to  arbitrary 
wing  plan  forma  undergoing  arbitrary  maneuvers.  The  boundary  values  for 
such  general  problems,  however,  usually  lead  to  the  development  of  double 
integral  equations  which  are  difficult  to  solve.  As  is  usual  in  such 
case’s,  there  are  many  special  types  of  plan  forma  and  maneuvers  which  lead 
to  boundary-value  problems  that  are  simpler  to  analyze.  An  important 
class  of  these  simplified  problems  is  that  arising  from  homogeneous  bound¬ 
ary  conditions. 

Let  cp(x,y,z,t)  be  a  solution  to  equation  (3).  In  certain  special 
cases  this  can  be  written  cp  =  (t)^cpQ(^,  in  which  case  cp  is 
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called  a  homogeneous  function  of  degree  n.  The  numher  of  variables 
affecting  9q  Is  only  three  as  compared  to  the  four  which  are  necessary 
to  determine  cp.  If,  therefore,  a  partial  differential  equation  can  he 
set  up  for  cpQ,  it  will  contain  one  less  mathematical  ** dimension"  than 
the  equation  for  cp  ,  Following  this  observation  it  is  necessary  to  pro¬ 
ceed  in  two  directions;  one  to  find  the  partial  differential  equation  for 
CpQ,  and  the  other  to  find  the  physical  problem  and  consequent  boundary 
values  leading  to  a  homogeneous  flow  field.  The  latter  path  will  be 
first  explored. 


First,  consider  an  example 
of  a  homogeneous  boundary— value 
problem.  Suppose  that  a  rectan¬ 
gular  flat  plate  at  an  angle  of 
attack  starts  suddenly  from  rest 
and  moves  forward  at  a  supersonic 
Mach  number  Mq.  At  "time”  ti 
the  initial  spherical  wave  gen¬ 
erated  by  the  forward  rlght-liand 
comer  has  traveled  outward  to  a 
radius  ti  and,  at  "time”  2ti, 
to  a  radius  2ti.  The  sketch 
Indicates  the  traces  of  these 
spheres  in  the  z  =  0  plane 
together  with  the  original  and 
present  position  of  the  wing 
leading  edge.  Let  the  points  Pi 
and  P2  be  located  on  the  same 
rays  through  the  origin  of  the 
circles  and  the  wing  comers.  The 
problem  is  to  find  the  pressures 
at  ?!  and  P2. 

It  is  apparent  that,  if 
every  dimens i oh  in  the  figure 
involving  P2  is  divided  by  2ti 
and  every  dimension  in  the  figure 
involving  Pi  is  divided  by  ti, 

the  two  figures  will  be  similar  in  every  respect  and  point  Pi  will  coin¬ 
cide  with  point  P2.  Since  the  vertical  velocity  Wq  is  constant  over 
the  plan  form,  a  simple  change  in  scale  has  made  the  boundary  conditions 
for  both  problems  identical.  But  this  means  that  the  solutions  at  Pi 
and  P2  are  identical  since  the  wave  equation  is  invariant  to  change  in 
scale.  Hence,  in  regions  of  a  rectangular  wing  unaffected  by  the  waves 
from  the  trailing  edge,  the  pressure  can  be  written 
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and  the  presaior©  Is  a  homogeneous  function  of  degree  zero.  A  generaliza¬ 
tion  of  this  example  is  contained  in  the  following  statement: 

(1)  The  pressure  in  any  region  affected  by  only  two  intersecting 
edges  of  a  straight-sided  flat  plate  traveling  at  a  uniform 
subsonic  or  supersonic  speed  is  homogeneous  and  of  degree 
zero  (i.e.,  satisfies  equation  (I6)), 

Consider  as  another  example  the  case  of  a  flat  rectangular  wing 
traveling  foirward  at  a  subsonic  or  supersonic  speed  and  rolling  about 
one  edge,  taken  to  be  coincident  with  the  x  axis.  The  argument  follows 
the  same  lines  as  before,  and  again  a  change  in  scale,  proportional  to 
the  time,  makes  ‘6h©  geometry  of  the  wing-^ave  combinations  identical  (in 
regions  affected  by  only  two  intersecting  edges)  for  different  times. 

The  boundary  values  over  the  wings  will  not  be  the  same,  however,  unless 
the  slope  of  the  Wq  distribution  is  adjusted  in  each  case.  But  Wq 
can  be  adjusted  by  reducing  it  an  amount  proportional  to  the  distance 
from  the  axis  of  rotation.  The  boundary-^alue  problems  are  then  similar 
for  different  values  of  time.  Finally,  therefore,  the  pressure  can  be 
written 


^  (i,y,z,t)  =7^(1,  |)  (IT) 

which  is  a  homogeneous  function  of  degree  one.  A  generalization  of  this 
example  is  expressed  as  follows: 

(2)  The  pressure  in  any  region  affected  by  only  two  intersecting 
edges  of  a  straight— sided  flat  plate  traveling  at  a  uniform 
subsonic  or  supersonic  speed  and  rotating  at  a  constant  rate 
of  pitch  or  roll  is  homogeneous  and  of  degree  one  (i.e., 
satisfies  equation  (I7)). 

It  should  be  noted  that  both  (l)  and  (2)  are  equally  true  for  the 
steady— state  case  when  all  transient  effects  have  disappeared.  In  super¬ 
sonic  wing  theory  they  lead  to  conical  and  quasi-conical  flows,  respec¬ 
tively,  while  in  the  subsonic  case  they  lead  to  flows  about  wings  having 
Infinite  chordwise  extent.  In  general,  homogeneous  flow  occurs  when  the 
bour).dary  conditions  after  a  change  in  scale  are  proportional  to  their 
original  values. 

Consider  next  the  modification  of  the  basic  partial  differential 
equation  (equation  (3))  under  the  assxmiption  that  the  flow  is  homogeneous. 
If  the  pressure  is  given  by  a  function  that  is  homogeneous  and  of  degree 
zero,  then,  by  equation  (6),  the  velocity— potential  function  will  be 
homogeneous  and  of  degree  one.  If  the  notation 
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^  ^  “  yo)  ^  -  zq 


(p  (x,y,2,t)  =  t«  (xq,  yo,  Zq) 


is  Tised,  then  equation  (3)  "becomes 


(l-^o^)®XoXo  +  (l-yo^)®yo7o  -  ^^0^0  *Xoyo 

^^0^0  ®x„z„  ~  ^Jo^o  ^y^z»  =  ° 


and  a  linear  partial  differential  equation  ■with  three  independent  vari¬ 
ables  is  therefore  obtained. 

In  the  general  theory  of  '  ^ 

partial  differential  equations  of  T  ° 

second  order ^  the  character  of  an  1 

equation  is  determined  from  the  _ . 

geometric  nature  of  a  related 
quadric  surface.  The  character  of 
equation  (19)  can  be  shown  from 
such  considerations  to  depend  on 
the  sign  of  the  expression  I-Xq^— 
yo^-Zo^-  It  ia  immediately 
apparent j  ho-wever^  that  within 
the  unit  sphere  in  the  x^^y^^ZQ 

space  the  sign  of  1— x:o^“-yo^-^o  ^  ' 

is  everywhere  positive  and  outside  .y 

the  sign  is  everywhere  negative.  t  ^ 

It  follows  that  outside  the  unit  I 

sphere  equation  (19)  is  hyperbolic 
and  inside  the  unit  sphere  it  is 

elliptic.  f  '^lA 

The  character  of  equation  (19)  Vy^  ^ 

is  of  particular  interest  since  the  A/ 

difficulties  inherent  in  the  deter—  - 

mination  of  the  solutions  can  be  ^  A 

estimated  without  actually  obtaining  A  Av 

the  solutions.  For  example,  con—  ^  ^ 

aider  the  two  configurations  shown 

in  the  sketch.  These  wings  started 

moving  at  t  =  0  with  the  foremost 

portion  of  their  leading  edges  on 

the  Jo  axis  and  have  by  now  traveled' forward  at  a  supersonic  speed  to 
attain  the  positions  represented  by  the  sketch,  the  unit  circle  being  in 
each  case  the  trace  of  the  primary  wave  from  the  vertex  on  the  z  =  0 
plane.  Outside  the  unit  sphere,  the  governing  equation  is  hyperbolic  and 
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the  hehavlor  of  the  flow  is  slinilar  to  that  In  steady-state  supersonic- 
wing  problems;  closed  solutions,  therefore,  can  he  expected  in  many  cases. 

Inside  the  unit  sphere,  on  the  other  hand,  the  character  of  equation  (19)  » 

is  elliptic.  The  canonical  form  of  an  elliptic  partial  differential  equa¬ 
tion  in  three  dimensions  is  Laplace’s  equation  which  is  the  governing 
equation  in  steady— state  suhsonic-^ing  theory.  Judging  from  the  com¬ 
plexity  of  subsonic,  three-dimensional  wing  problems,  it  can  be  concluded 
that  a  closed  solution  within  the  unit  sphere  will  be  very  difficult  to 
find. 


It  is  Instructive  to  notice  that  this  entire  development  has  a  direct 
analogue  in  the  study  of  three-dimensional ^  steady  state,  supersonic  wings. 
In  that  case  the  original  equation  is  the  three-dimensional  wave  equation 

'Pxx-'Pyy-'Pzz=  0  (20) 

By  considering  the  velocity  potential  to  be  homogeneous  and  of  degree  one, 
Busemann  in  reference  T  was  able  to  introduce  the  transformations 


7o  = 


y 

X 


(p(x,y,z)  =  X  $  (yo,Zo) 

and  transform  equation  (20)  to  the  form 


(1-yo^)  Voyo  +  (1-Zo^)  ^Zq^o  "^^0^0  ^yo^o  =  ° 

which  is  the  two-dimensional  form  of  equation  (l9).  Flows  governed  by 
equation  (21)  have  become  known  as  conical  flows.  A  study  of  equation  (21) 
shows  it  to  be  elliptic  inside  and  hyperbolic  outside  the  mit  circle. 

In  this  case,  however,  the  equation  has  only  two  independent  variables  so 
that  once  the  equation  has  been  transformed  to  the  two-dimensional  form 
of  Laplace’s  equation®  solutions  are  not  difficult  to  find. 


®The  Tschapligin  transformation  transforms  equation  (21)  into  Laplace’s 
equation. 
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Boundary— Value  Problems  Involving  Nonlnterac ting  Surfaces 


Tbe  simplification  of  the  four-dimensional  theory  brought  out  by 
the  introduction  of  homogeneous  flow  was  more  apparent  than  real  since  ' 
the  resulting  partial  differential  equation,  although  containing  one 
less  dimension,  was  unwieldy.  Another  class  of  wing  problem  which  is 
simplified  both  in  theory  and  in  practice  by  reasoning  from  physical 
knowledge  of  the  flow  behavior  is  that  in  which  the  wing  has  a  super¬ 
sonic  edge  (i.e.,  an  edge  which  la  traveling  with  a  supersonic  normal 
component  of  velocity) . 

When  the  acoustic  plan  form  is  affected  only  by  a  supersonic  edge, 
it  is  obvious  that  the  flow  on  the  upper  surface  of  the  wing  is  independ¬ 
ent  of  that  on  the  lower  surface.  Hence  the  solution  to  such  problems 
can  always  be  written  in  terms  only  of  sources  as  follows: 


dS 


(22) 


where  Scp/Szi  =  Wu(xi,yi)  = 

\  ^u  being  the 

local  slope  of  the  surface  in 
the  direction  of  Vq.  Since 
equation  (22)  is  equally  valid 
for  symmetrical  nonlifting 
surfaces  and  lifting  plates, 
its  value  and  simplicity  is 
evident. 

If  the  wing  plan  form  is  ' 
further  specialized  hy  having 
not  only  supersonic  leading 
edges,  hut  also  having  a 
straight  trailing  edge  perpen¬ 
dicular  to  the  direction  of 
motion,  additional  simplifica¬ 
tions  can  he  used,*^  Consider, 
for  example,  the  two-dimensional 
wing  (a)  in  the  sketch.  Let 
this  wing  have  an  angle  of 
attack  a(t)  which  varies  with 
time  in  an  arbitrary  manner. 

There  results  from  such  an 
angle-of -attack  variation  a 
certain  lift  which  also  varies 
with  time.  Hence,  if  L*  represents 


the  total  lift  on  an  airfoil  of 


*^The  following  method  simply  extends,  to  include  the  effects  of  unsteady 
motion,  a  theorem  given  hy  Lagerstrom  and  Van  Dyke,  ,  (See  reference  2.) 
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very  high  aspect  ratio  and 
slope,  then® 


represents  the  section  lift-curve 


—  =  c  ciJ  (c,t)(span) 
qa 


Next,  it  is  clear  by  reason  of  symmetry  that  the  total  lift  on  vings  (h) 
and  (c)  in  the  sketch  are  equal.  Then,  since  the  analysis  is  based  on  a 
linear  partial  differential  equation,  by  superposition  principles  the 
total  lift  on  vlng  (b)  or  (c)  is  Just  half  of  that  on  wing  (a).  In 
another  sense,  the  lift  coefficient  for  the  whole  wing  based  on  the 
deflected  area  is  the  same  for  all  three  cases.  A  suitable  superposition 
of  wings  (a),  (b),  and  (c)  will  give  wing  (d),  which  then  has  the  same 
lift  coefficient  based  on  deflected  area.  Finally,  because  of  the  super¬ 
sonic  stream,  wing  (e)  can  be  obtained  from  (d),  hence  it  also  has  the 
lift  coefficient  common  to  the  other  wings.  It  is  of  course  necessary 
that  the  variation  of  a  with  time  be  the  same  in  each  case. 


The  preceding  process  can  be  extended  one  step  farther  to  the 
development  of  the  lift  due  to  a  single  deflected  element.  By  consider¬ 
ing  the  sketch  it  can  be  seen  that 


®The  primes  on  quantities  indicate  that  two-dimensional  values  are  taken, 
or,  that  a  high-aspect-ratio  wing  is  considered  and  tip  effects  are 
neglected.  _ 
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Finally,  if  a  wing  is  composed  of  a  distribution  of  these  elements,  then, 
for  an  axial  system  centered  at  the  apex  of  the  wing  leading  edge,  there 
results 

Jg  J  L  J 

where  cq  is  the  maximum  chord  (aee  sketch).  In  the  development  of  equa¬ 
tion  (23)  each  element  is  asaiamed 
to  have  the  same  variation  of 
motion  with  time. 

Notice  that  when  the  wing  is 
a  flat  plate  flying  at  a  steady 
speed  so  that  all  transient  effects 
have  disappeared,  c^  *  is  independ¬ 
ent  of  t  and  of  the  chord  length, 
being,  in  fact,  equal  to  U/p.  Then 
equation  (23)  becomes 

Cl  =  ^//a(x,y)d^dy  .  f 
s 

where  a  is  the  average  angle  of  attack  of  the  wing.  This  result  has 
already  been  obtained  in  reference  2.  When  a(x,y)  is  independent  of 
X  (as  for  a  flat  wing  sinking  or  rolling),  equation  (23)  becomes 

b 

/I  ^^(y)  o  (c,t)  dy  (2h) 

b  a 

where  c  is  the  local  chord  which  is,  in  general,  a  fixaction  of  y. 
Equation  (24)  simply  indicates  that  longitudinal— strip  theory  is  exact 
for  such  wings. 

Finally,  notice  that  the  calculation  of  the  uasteady  lift  on  three- 
dimensional  wings  with  supersonic  leading  edges  and  straight  trailing 
edges  perpendicular  to  the  free— stream  direction  has  been  reduced  to  an 
Integration  involving  the  relatively  simple  results  for  a  two-dimensional 
wing  undergoing  the  same  unsteady  motion.  For  example,  the  lift  on  a 
flat  triangular  wing  rising  and  sinking  with  a  harmonic  motion  can  be 
computed  from  a  single  integration  of  the  results  presented  in  refer¬ 
ence  3.  Such  a  calculation  can  be  carried  out  numerically  quite  rapidly. 

Two-Dimensional  Boundary-^Talue  Problems 

The  simplification  brought  about  when  the  flow  is  independent  of 
one  dimension  is  again  obvious.  In  such  cases,  the  three-dimensional  wave 
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equation  (3)  reduces  immediately  to  the  two-dimensional  wave  equation. 
Typical  examples  of  this  type  of  problem  can  be  constructed  by  consider¬ 
ing  flat  plates  which  start  suddenly  at  t  =  0  and  travel  thereafter  at 
constant  supersonic  velocities.  Two  examples,  one  a  comer  of  a  rectan— 
giilar  wing  and  the  other  a  triangular  flat— plate  wing  are  shown  in  the 
sketch.  After  time  t  =  0,  the  edges  of  the  wings  send  out  cylindrical 
waves  and  the  outer  boundaries  of  these  waves  at  time  t  are  shown  as 
dashed  lines  parallel  to  the  edges  in  question. 


Leading  edge  now 


Mach 

waves 


Leading  edge 
at  t=0 


Since  points  in  regions  1  and  2  are  affected  only  by  a  single  edge,  the 
wave  phenomena  in  these  regions  are  cylindrical,  and  the  physical  quan¬ 
tities  are  in  both  cases  Independent  of  distance  parallel  to  the  edge 
which  acts  as  their  generator.  Hence,  the  flow  field  in  these  regions 
may  be  regarded  as  two  dimensional.  (Region  3^  incidentally,  is  independ¬ 
ent  of  distance  in  both  the  x  and  the  y  directions  and  is,  therefore, 
one-dimensional. ) 

Solutions  to  the  two-dimensional  unsteady  problems  are  sometimes 
especially  easy  to  find  because  of  the  analogy  they  have  with  three- 
dimensional,  steady-state,  liftlng-eurf ace  problems,  (See,  e.g,,  refer¬ 
ence  8.)  For  example,  consider  an  Infinitely  long  unyawed  wing  which 
starts  from  rest  and  travels  forward  at  a  velocity  Y  which  may  or  may 
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not  te  a  function  of  time.  The  trace  of 
this  wing  in  the  Xyt  plane  is  like  that 
shown  in  the  sketch,  (in  the  sketch  shown^ 
the  wing  velocity  is  always  less  than  the 
speed  of  sound  and  varying.)  The  "boundary 
conditions  are  that  is  specified  over 

the  shaded  area  and  the  loading  Acp^  is 
zero  everywhere  except  within  the  shaded 
area.  But  if  x  is  replaced  by  y  and 
t  by  Xj  these  boundary  conditions  are 
exactly  the  same  as  those  for  a  plate  of 
known  camber  and  angle  of  attack,  with  a 
plan  form  as  indicated  by  the  shaded  area, 
placed  in  a  free  stream  directed  along  the 
positive  X  axis  at  a  Mach  number  equal 
to  y  2.  The  solution  for  the  one  problem 
may  be  used,  therefore,  as  a  solution  to 
the  other  with  only  a  change  in  notation. 


Traces  of 


characteristic  cones 


Boundary  Conditions  for  Very  Slender  Wings 


When  the  wing  plan  form  is  slender  in  the  sense  that  its  length  in 
the  streamwise  direction  is  large  compared  to  its  span,  an  estimation  of 
the  loading  on  it  can  be 
obtained  by  neglecting  in 
the  partial  differential  equa¬ 
tion  the  gradient  of  the  Induced 
velocity  component  in  the  stream 
direction.  Thus,  if  the  wing 
is  moving  in  the  negative  x 
direction,  equation  (3)  reduces 
to 


A - 


^tt 


(25) 


which  la  again  the  wave  equa¬ 
tion  but  in  two  space  dimensions. 
Since  equation  (25)  is  independent 
of  X,  study  can  be  made  independ¬ 
ently  of  the  boundary  conditions 
on  each  spanwise  strip.  This  is 
an  extension  of  steady-state 
slender-wing  theory.  (See,  e.g., 
reference  9.)  The  sketch  shows  a 
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tn>ical  section  in  the  yt  plane.  If  the  wing  is  a  flat  plate  starting 
from  rest  and  traveling  forward  at  a  ■uniform  speed,  the  boundary  condi 
tion  is  that  (Pz  should  he  constant  over  the  shaded  area  in  the  sketch 
and  Acp+  should  he  zero  everywhere  except  across  this  area.  The  solution 
to  this  problem  leads  directly  to  the  solution  for  the  thin  triangular 
wing.® 


Two-Dimensional  Ikisteady  Incompressible  Flow 


The  analogy  between  two-dimensional  unsteady  and  three-dimensional 
steady  flow  provides  an  Interesting  viewpoint  for  two-dimensional  unsteady 
problems  in  incompressible  flow.  In  this  case  the  boundary  conditions 
are  independent  of  the  lateral  coordinate  and  equation  ( 1 )  becomes 

1 

'Pxx+'Pzz  =  ^  'Pt't* 

Since  the  speed  of  pressure  propagation  ao  la  assumed  infinite,  the 
■basic  equation  thus  'becomes 


Leading 

edge 


Trailing 
edge 


Vortex 

sheet 


vhere,  however,  time  still  appears 
in  the  boundary  conditions  and  in 
the  expression 


q 


at* 


which  is  eqvil valent  to  equation  (6). 
The  technique  involved  in  the  solu¬ 
tion  of  such  problems  is  analogous 
to  that  in  which  one  term  of  the 
three-dimensional  steady-state 
eqviation  vanishes  by  virtue  of  the 
fact  that  free-stream  Mach  number 
approaches  one. 


If  the  two-dimensional  wing 
starts  from  rest  and  travels  for¬ 
ward  at  a  speed  V,  the  trace  of 
the  wing  is  as  shown  in  the  sketch 
between  this  problem  and  the  more  general  case 


The  essential  difference  - -  -  -  ^ 

of  two-dimensional  compressible  flow  lies  in  the  fact  that  the  traces  of 

®It  is  obvious  that  the  analogous  problem  in  steady— state  wing  theory  is 
that  of  a  low-aspect— ratio,  rectangular,  flat  plate  in  a  froe  stream 
having  a  Mach  number  equal  to  _ 
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the  characteristic  cones  are  normal  to  the  t*  axis.  The  ■boundary- 
conditions  are  therefore  satisfied  along  lateral  strips  and^  in  lifting- 
surface  terminology^  the  analysis  corresponds  to  Slender-^ing  theory. 

These  latter  methods  are  toII  estahlished  and  in  reference  9,  for  example, 
the  manner  in  which  the  Kutta  condition  is  imposed  is  discussed  in  some 
detail.  The  trailing  vortex  sheet  for  the  lifting  wing  has  the  same 
distribution  of  vorticlty  that  exists  behind  the  -unsteady  two-dimensional 
airfoil  and  the  rolling  up  of  the  vortex  sheet  can  be  studied  from  either 
standpoint .  ' 


APPH  CATIONS 


Starting  Lift  of  a  Wing 


One  of  the  simplest  and  yet  most  general  results  which  can  be  derived 
on  the  basis  of  the  present  theory  is  the  initial  value  of  pressure  on  a 
wing  surface  starting  impulsively  from  rest  with  a  velocity  Vq* 
discussion  will  be  made  for  a  wing  without  thickness  although  the  results 
will  be  seen  to  apply  to  the  thickness  case  as  well.  •; 


Consider  a  surface  with  a 
plan  from  as  indicated  in  the 
sketch.  The  acoustic  plan  form 
of  a  point  P(x,y)  on  the  sur¬ 
face  is  a  small  circle  of  radius 
t.  Since  no  point  on  the  wing 
outside  this  circle  can  influence 
the  pressure  at  P,  the  upper 
surface  is  independent  of  the 
lower  surface,  except  for  a  band 
of  width  t  around  the  edge  of 
the  wing.  It  is,  therefore, 
evident  that  the  boundary— value 
problem  to  be  solved  has  been 
treated  in  the  section  Boundary- 
Value  Problems  Involving  Noninter— 
acting  Surfaces,  The  solution 
follows  directly  from  equation  (22) 
and  can  be  written 
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Using  a  polar  coordinate  system  defined  by 


X  —  =  r  cos  9 
y  -  yi  =  r  sin  9 
dxi  dyi  =  rdrde 


there  results 


so  that 


^  SJ 


Wu(i,y) 

2jt 


dr  -  -  twu(x,y) 


^  Wu(x,y) 

V  q  ^oMoV'Sty^^Q  VoMo 


If  a  =  -  Wu(x,y)/V'o  is  the  local  slope  of  the  wing,  the  expression  for 
load  coefficient  becomes 

f _  Mx,y)  /ogx 

Vq  4o  ■  “o 


The  starting  value  of  lift  coefficient  can,  therefore,  he  written 


(27) 


where  a  la  the  average  angle  of  attack  of  the  surface  defined  hy 

I-iy’^adxdJ 

S  being  the  area  of  the  wing  plan  form. 


^°The  mean  value  theorem  gives 

CD  =  rr  — 

2jt  JJ  Tq 

where  |  and  tj  lie  somewhere  in  3a.  HencSj  as  t  approaches  0, 
I  and  T]  approach  x  and  y,  respectively. 
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Supersonic  Steady^-State  Lift 


In  order  to  illustrate  some  of  the  TDasic  concepts^  consider  next 
an  example  for  vhicli  the  result  is  well  known.  The  problem  will  be  to 
derive,  by  the  methods  outlined  in  this  report,  the  steady-state  loading 
.on  a  two-dimensional  flat  plate  traveling  at  a  constant  supersonic  speed. 

7 

Since  the  upper  and  lower  surfaces  are  noninteracting,  the  solution 
stems  from  equation  (22).  In  the  plane  of  the  wing  this  equation  becomes 


Leading  edge  now 


The  accompanying  sketch  shows  the  positions  of  the  wing  in  the  xy 
plane.  The  wing  has  constant  speed  for  t>0.  The  point  P(x,y)  is 
chosen  on  the  wing  and  ahead  of  the  wave  which  started  at  time  zero; 
therefore,  P(x,y)  lies  in  the  region  which  has  attained  its  steady-state 
value.  Further,  the  value  of  is  constant  over  the  acoustic  plan 

form  in  such  a  region.  This  constancy  reduces  the  problem  to  one  of 
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integrating  [ (x— xi)^+(y-yi)^]  over  the  ellipse  representing  the 

acoustic  plan  form. 

The  eqmtlon  for  can  be  determined  from  equations  (l4)  and 

(15).  In  this  cassj  equation  (15)  becomes  simply 

Xi  =  -  MqT 

Eliminate  t  between  this  equation  and  equation  (l4)  and  there  results 
in  the  z  =  0  plane 

(x-xi)^  +  (y-yi)^  =  ^  +  (29) 

That  equation  (29)  is  the  equation  of  an  ellipse  with  one  focal  point 
at  XjY  can  he  readily  verified.  It  is  more  convenient^  however^  to 
change  to  a  polar  coordinate  system  with  origin  at  P.  Hence  set 


X— Xi  =  r  cos  B 
y— yi  =  r  sin  B 
dyi  dXi  =  rdrd6 


Then  equation  (29)  becomes 


Mgr  =  Mot  X— r  cos  B 


or 


r  = 


X  +  Mot 

Mq  +  cos  B 


and,  therefore,  equation  (28)  becomes  simply 

X  +  Mot 


,p  . f  ief ^ 


Wu  /  ^  .  V  de 

=  -  ^  {x+  Mot)  /  - - 

Mn  +  < 


Mo  +  COS  B 


The  integral  is  not  difficult  to  evaluate  so 


<P  =  -  —  (x+M  t)  ^ 
2it  °  ' 


(30) 
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and,  finally,  by  equation  (6) 

Ap  _  4  /  WuMo  \  it-a  /  , 

q  "  voMo  y  =  (31) 

which  is  the  familiar  Acker et  value  for  the  loading  on  a  two-dimensional 
flat  plate.  The  lift  coefficient^  of  course_,  follows  immediately  as 


4a 

VMo^-I 


(32) 


Stability  Derivative  C7 

^r 


Consider  the  case  of  a  hlgh-aspect-ratio  (so  that  the  tip  effects 
can  he  neglected)  rectangular  wing  yawing  at  a  uniform  rate  r  radians 
per  second  and  traveling  forward  along  its  flight  path  at  a  supersonic 
velocity  Vq.  The  eqiiatlon  of  the  dealing  edge  of  the  wing,  using  the 
coordinates  indicated  hy  the  sketch,  is  given  hy  the  relation 


X 


and  since  ilr  =  rt*  this 
becomes 

^  =(7  -  tan  rt*  (33) 

The  boundary  conditions 
(when  the  stability  derivative 
is  to  be  referred  to  the  body 
axis)  are  that  w^,  the  ver¬ 
tical  induced  velocity  in  the 
region  occupied  at  any  given 
time  by  the  plan  form,  shall  be 
constant.  Further,  the  steady- 
state  condition  is  determined 
when  the  Intersection  of  the 
leading-edge  trace  with  the 
inverse  sound  wave  forma  a 


Wing 


closed  curve.  Honce,  the 

acoustic  plan  form  is  given  by  eliminating  t 


from  the  equations 
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(x-xi)2+  (y-yi)2=  (t-T)2 


(’■-» 


tan  rr* 


(31^) 


It  is  usual  to  assume  r  to  be  small  and  to  retain  only  the  first 
power  of  r  in  a  aeries  expansion.  Therefore,  since  it  can  be  shown 
that  the  equation  of  the  acoustic  plan  form  is  given  correctly  to  the 
first  power  in  r  by  assuming  rt’  =  tan  rt*,  equation  (3^)  becomes 


where,  as  in  the  preceding  example,  the  polar  coordinates 

X— Xi  =  R  cos  9 
y— 7i  =  H  sin  9 
dxi  dyi  =  EdRd0 

have  been  introduced  together  with  the  notation 

r*  =  — 

So 

solving  equations  (35)  ^or  R,  there  results 


E  = 


-A^  -  4  Bpr*  sine 

2r*  sin  9 


(36) 


(37) 


where 


Aq  =  cos  e  —  +  Mq  —  sin  6 

Bo  =  rny-x-Mot 

and  expanding  in  powers  of  r*  gives,  neglecting  powers  higher  than  the 
first. 


or 


R  = 


x+Mot 


Mo+  cos  9 


ty 


|_Mo+  cos  9 


(x+Mot)(y+t  sin  9) 

(Mo+  COS  9)^ 


(x+Mot) 


(Mo+  cos  0)3 


sin  9 


31 


MCA  TN  2256 


The  solution  for  the  velocitj  potential  is  given,  just  as  in  the 
preceding  example,  hy  the  equation 

cp  =  _  /  ecL0 

and  since  , 


o 


sin  Q 


(Mq  +  cos  0)® 


sin  0 


(Mq  +  cos  0) 


d0  =  0 


'q  Mq  +  cos  0  v^MqS-I 

n  2Jt 

/  d0  ^  2rtMo 

(Mo  +  cos  0)2  (Mo^-l)  3/2 


then  finally 


cp  - - y-  ■  (Mpt  +  X  +  r^  --g--.. 

-v/Mo2-1\^  Mo^-l 


The  loading  can  he  obtained  from  equation  (38)  hy  using  equation  (6) 


hence 


Ap  4a  4a 

—  =  —  +  - 

<3  J  M^2_i  (Mo2-1) 


\^oJ 


This  result  has  "been  obtained  previously  in  reference  10  by  a  different 
method. 


Moving  Two-Dimensional  Vortex 


Moving 

vortex 


A  further  interesting  application  of  the  m/n  n  t  Starting 

methods  presented  here  is  the  development  of  ^  vortex 

the  moving  rectilinear  vortex,  where  a  vortex 

is  represented  by  a  line  of  constant  clrcula—  T  T  x 

tion.  The  moving  vortex  is  associated  with  a  ^  ^o*^  ^ 

lifting  surface^^  which  starts  Impulsively  at 
t  =  0  and  travels  at  a  imiform  velocity  away  • 
from  a  fixed  axial  system.  Since  a  vortex 

cannot  end  in  a  fluid  (Heloiholtz 's  fimdamental  result),  the  moving  vortex 
must  be  adjoined  by  two  trailing  vortices  and  they,  in  turn,  by  a  final 
starting  vortex  (see  sketch)  which  remains  in  the  vicinity  of  the  position 


Associated  in  the  sense  that  a  distribution  of  these  vortices  will 
satisfy  the  boimdary  conditions  for  the  lifting  siarface.  The  moving 
vortex  is  I 


mm 
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at  which  the  motion  commenced.  For  the  present^  the  assumption  is  made 
that  the  trailing  vortices  are  far  enough  away  so  that  they  do  not  affect 
the  flow  around  the  center  of  the  moving  and  starting  vortices  (i.e.j  the 
flow  is  two  dimensional);  and  second,  that  the  starting  vortex  remains 
fixed  on  the  line  along  which  the  moving  vortex  first  appeared.  With 
these  assumptions,  the  boundary  conditions  can  be  determined  from  the 
definition  that  the  vortex  system  described  is  equivalent  to  a  distribu¬ 
tion  of  doublets  of  constant  strength  over  the  region  of  the  z  =  0 
plane  bounded  by  the  moving  vortex,  the  stationary  vortex,  and  the  vor¬ 
tices  Joining  them  at  their  extremities.  Such  a  sheet  of  doublets  gives 
rise  to  a  constant  potential  difference,  Acp,  across  the  sheet,  while  the 
normal  derivative  is  continuous  so  that  A(S'P/Sz)  =  0.  The  solution  is 
then  given  by  equation  ( I3 ) .  Thus 

cp  =  cpi  +  £p2  (*^0) 

where 


and 


1  rr  i 

J c  \bt  /  bz 


Til©  evaluation  of  cp^  imposes 
no  difficulties  except  in  integra¬ 
tion  but  the  calculation  of  9 2 
requires  some  special  consideration. 
Sine©  A9  is  a  constant, 
is  zero  everywhere  except  along 
the  line  of  the  moving  vortex  and 
there  it  is  infinite.  In  order  to 
evaluate  such  an  indeterminate 
form,  assume  that  the  potential 
difference  rises  linearly  to  its 
value  Acp  in  the  time  interval  g 
(see  sketch).  Then  the  value  of 
^Acp/St  is  AP/g  in  the  narrow 
strip  between  the  two  ellipses 
formed  by  the  intersection  of  the 
inverse  sound  wav© 

(x-xj^  +  (y-yi)^  +  z^  =  (t-T)^ 
with  the  lines 

Xl 


=  -MoT 

=  -  Mo(t-  e) 
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■When  the  equations  for  these  two  ellipses  are  formulated  and  placed  in 
equation  (42)  there  resiilts  an  equation  which  can  he  represented  by  the 
expression 


'P2 


_  A(p 

hire 


F(t)-^’(t-e) 


where 


F(t)  = 


q  ^0 


and  So  ia  the  aam©  acoustic  plan  fom  used  in  equation  (4l).  But  since 


11m 
e  — >0 


F(t)-^(t-e)  ^  ^ 
e  bt 


the  equation  for  ^2  "becomes 

hjt  btJ  Ja  bz 

^a 

The  calculation  of  cpi  and  cp2  over  the  acoustic  plan  form  given 


AT  S  r  r  1  Sro 


hy  the  area  within  the  ciurve 

(x-xi)^  +  (y-yi)^  +  z2  =  ^ 

or  between  this  curve  and  the  yi 
axis  gives  for  the  regions  defined 
by  the  sketch 


(43) 


Region  (l) 


CD  =  ^ 
^  2jt 


arc  tan 


(x+Mot)  </t^— x^— Moving 


z(MQX+t) 


vortex 


•  arc  tan  ^ 


Region  (2) 


2  2 
-x  ~z 


(44) 


<p  = 

^  2}f 


arc  tan 


Region  (3) 


3*^ 
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It  la  to  be  noted  that  eqixatlon  (hh)  gives  the  potential  everywhere 
within  the  cylinder  In  which  disturbances  can  be  affected  in  subsonic 
flow  as  well  as  over  a  large  portion  of  the  disturbed  region  in  super¬ 
sonic  flow. 

It  is  apparent  that  the  results  given  in  equations  (44),  (45),  and 
(46)  can  also  be  obtained  independently  by  the  methods  indicated  in  the 
section  entitled  Two-Dimensional  Boundary— Value  Problems.  From  this 
viewpoint,  the  problem  of  the  unsteady  two-dimensional  vortei  pair  has 
for  its  analogue  in  lifting-surface  theory  the  case  of  a  triangular  wing 
yawed  so  that  one  side  is  parallel  to  the  free-etream  direction  and  the 
other  side  is  either  a  subsonic  or  a  supersonic  edge,  depending  on  the 
case  being  considered.  Since  Acp  is  specified  a  constant  over  the 
entire  plan  form,  the  solution  is  relatively  easy.  In  reference  11,  for 
example,  the  solution  for  the  subsonic  case  is  given  in  terms  of  the 
variables  used  in  steady-state  analysis. 

The  foregoing  results  for  Mo<l  reduce  to  more  familiar  expres¬ 
sions  when  their  asymptotic  expressions  for  large  values  of  time  are 
determined.  A  study  of  the  two  terms  in  equation  (44)  shows  that  as  t 
becomes  large  the  expression  becomes 

^Ao/jtz  .  xN 

2^  (,2F| 

This  result  provides  the  velocity  potential  for  a  stationary  vortex  in 
an  incompressihle  medium,  the  jump  in  the  value  of  the  potential  occurr¬ 
ing  along  the  plane  of  the  vortex  motion  (see  sketch). 

If  the  axes  are  fixed  in 
the  moving  vortex,  nev  coo2>- 
dinates  I,  t  need  to  he 
introduced  as  defined  hy  the 
equations 

g  =  X  +  Mot 

^  =  z 

T  =  t 

As  T  becomes  large,  equa¬ 
tion  (^4*4)  reduces  to  the  form 

'P  S’  iTl  t  h  -•==) 

2  1^1  C-/L-Mo2  / 

This  result  agrees  with  that  of  G-lauert  (reference  12)  for  a  vortex 
obtained  from  the  steady-state  linearized  theory  of  compressible  flow, 
so  that  in  the  vicinity  of  the  moving  vortex  the  induced  velocity  field 
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is  identical  to  the  perturbation  field  of  a  fixed  vortex  in  a  subsonic 
free  stream. 


It  is  interesting  to  inspect  the  vertical  velocitj  induced  by  the 
vortex  pair  in  the  z  =  0  plane.  From  equations  (4h)  and  (45)  it 
follows  that  for  x^<t^  the  vertical  velocity  w^  for  both  subsonic 
and  supersonic  speeds  is  given  by  the  relation 


Acpr/WoX+t  t\  1  1 

mq 

2jt  x(,x+MQt;  ^ 


(1^7) 


The  sketch  shows  the 
variation  of  w-^;^  along  the 
X  axis  for  a  moving  vortex 
traveling  at  both  subsonic 
and  supersonic  speeds ,  For 
both  cases  in  the  vicinity 
of  the  starting  vortex  w^ 
behaves  as  it  does  near  a 
two~4imensi  onal,  s teady-nstate , 
fixed  vortex,  having  large 
values  of  downwash  on  one 
side  and  large  values  of 
upwash  on  the  other.  The 
velocity  distribution  around 
the  moving  vortex  is  quite 
different,  however,  depend^ 
ing  on  whether  its  Jyfech 
number  is  greater  or  less 
than  one.  In  the  subsonic 
case  the  distribution  is 
similar  to  that  for  the 
starting  vortex,  but  in  the 
supersonic  case  the  velocity 
is  just  a  pulse  existing  at 
the  position  of  the  vortex. 
Again  in  both  cases  the 
velocity  falls  to  zero  at 
the.  primary  wave  front  (i.e., 
at  X  =  ±  t  in  the  sketch) . 


X 


X 


A  final  remark  can  be  made  concerning  equation  (kk)  in  connection 
with  incompressible  theory.  In  this  case,  the  value  of  Uq  becomes 
Infinitely  large  and,  if  equation  (ii-4)  is  written  in  the  form 
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Acp  .  (x+Vot*)  vao^'t'^3c2-2  2  x 

arc  tan  - - - 2C_SJ - arc  tan  - -  - - — - 

^r^-ov) 


it  is  apparent  that^  as  aQ->oo  the  velocity  potential  ■becomes 


Cp  ^  ^  - 

2jt  V  z 


x-hVot* 


—  arc  tan 


Equation  (48)  gives  the  velocity  potential  for  a  moving  and  a  fixed 
vortex  pair  in  an  Incompressihle  fluid. 


Moving  Three-Dimensional  Vortex 


In  three-dimensional  subsonic ^2  unsteady  lift  theory  an  Important 
role  is  played  by  the  loop  vortex  composed  of  a  moving  bound  vortex  of 
fixed  span^  a  stationary  starting  vortex  of  opposing  strength^  and  two 

trailing  vortices  connecting  the 
end— points  of  the  other  two  and 
completing  the  vortex  loop.  The 
boundary  condition  for  such  a 
configuration  requires  that  the 
jump  in  potential  Acp  la  a  con¬ 
stant  over  the  interior  of  the 
loop.  The  calculation  of  the 
Induced  field  can  be  resolved 
essentially  into  the  determina¬ 
tion  of  the  velocity  potential 
for  two  semi— infinite  vortices 
and  a  trailing  vortex  (see 
sketch)  aince^  by  superposition 
methods,  two  of  these  latter 
configurations  can  be  combined 
with  the  two-dimensional  figure 
of  the  last  section  to  give  the 
required  loop.  The  present 
example  is  concerned  with  the  calculation  of  the  velocity  potential  in 
the  tip  region  of  the  semi— inf ini te  case. 

As  in  the  two-dimensional  case  just  discussed  the  potential  can  be 
divided  into  two  parts  and  T  2  such  that 


12 


The  extension  to  Include  the  three-dimensional  supersonic  loop  is 
obvious. 
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9  =  9i  +  «P2 


(p  Aq>z  rr  dxidyi 

"  hitJJ  [(x-Xi)2+(y-y,)2+z2]3/2 

92  =  +  ^  —  A  r 

Im  j  Jg  [(x-Xi)^+(y-7i)^+z^] 

d 


where  Sg^  is  the  acoustic  plan  form,  "being  the  area  bounded  by  the 
lines 

(x-xi)^  +  (y-yi)®  +  z®  = 

71  =  0 
Xi  =  0 

It  is  convenient  to  use  the  notation 


k  =  x+Mot' 

I  =  3mci 

^  =  7-71 

Then  the  potential  can  he  written 


cp  =  r  ^ _ 

-LV^t2_x2-z2 


zAcp_^  _ 

UT”  St  J_ 


d,/ 

X 


^  (k-Mov4:^+P^z^+P^n^) 


dl 


(I^+T1^Z^)=/2 


:  2 


dTi 


/ 


^  (k-Moyt^+p^z^+p^Ti^) 


iW+^ 
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which  becomes,  after  one  integration. 


f  ^  1  a,  * 

\T1^+Z^  k^+T)^+Z^/  J 


zA<p 

TT. 


^t2-x2-z2 


Mq  ^  _ Mo^ _  1 

k^+P+ip  (k^+z^+T)^)  V'ks+pszS+pTjs 


Combining  and  integrating  again  there  is  obtained,  for  Mp<l, 


cp  _  Acp  r 

4)t  L 


arc  tan 


(x+MQt)  //t^— 
z(MoX+t) 


—  arc  tan 


xJ 


arc  tan 


y(x+Mot) 


(^^9) 


which  is  the  potential  within 
the  primary  spherical  wave 
shown  in  the  sketch,  Notice 
that  across  the  shaded  area  the 
jump  in  potential  given  by 
equation  (l^-9)  is  equal  to  A9 
and  elsewhere  the  potential  is 
continuous.  Along  the  portion 
of  the  sphere  t^  =  y^'+z^+x^ 
for  which  y>0,  the  potential 
is  equal  to  that  derived  for 
the  two-dimensional  vortex;  and 
on  the  rest  of  the  sphere, 

<p  equals  0* 


Ames  Aeronautical  Laboratory, 

National  Advisory  Committee  for  Aeronautics, 
Moffett  Field,  Calif.,  October  12,  I95O. 
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